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The steady boundary layer flow and heat transfer of a nanofluid past a nonlinearly permeable stretching/ 
shrinking sheet is numerically studied. The governing partial differential equations are reduced into a 
system of ordinary differential equations using a similarity transformation, which are then solved 
numerically using a shooting method. The local Nusselt number and the local Sherwood number and some 
samples of velocity, temperature and nanoparticle concentration profiles are graphically presented and 
discussed. Effects of the suction parameter, thermophoresis parameter, Brownian motion parameter and 
the stretching/shrinking parameter on the flow, concentration and heat transfer characteristics are 
thoroughly investigated. Dual solutions are found to exist in a certain range of the stretching/shrinking 
parameter for both shrinking and stretching cases. Results indicate that suction widens the range of the 
stretching/ shrinking parameter for which the solution exists. 

The study of boundary layer flow and heat transfer over a stretching sheet has gained vast interest among 
researchers due to its various applications in industrial and engineering processes for example in manufac- 
ture and extraction of polymer and rubber sheets. Since then, various aspects of stretching sheet problems 
have been investigated by several researchers (see Kumaran et al. 1 , Hayat et al. 2 , Ishak et al. 3 , Mahapatra et al. 4 , 
Fang et al. 5 , Makinde and Aziz 6 , Rana and Bhargava 7 , Mahmoud and Megahed 8 , Rahman and Eltayeb 9 , Khan and 
Shahzad 10 , Mukhopadhyay 11 etc.). It is well known that Khan and Pop 12 presented the first paper on stretching 
sheet in a nanofluid. 

Since the past few years, much interest was focused on convective heat transfer in a nanofluid. The base fluids 
such as water, oil and ethylene glycol used in many industrial processes such as in power generation, chemical 
processes and heating or cooling processes are poor heat transfer fluid due to its poor heat transfer properties with 
low thermal conductivity. One can improve this by suspending the solid nanoparticle into the base fluid in order 
to increase the thermal conductivity. According to Daungthongsuk and Wongwises 13 , the poor heat transfer 
properties of the base fluids was identified as a major obstacle to the high compactness and effectiveness of heat 
exchangers. The essential initiative is to seek the solid particles having thermal conductivities several hundreds of 
times higher than those of base or conventional fluids. It seems that Choi 14 was the first to call the mixture of the 
base fluids with the solid nanoparticle as nanofluid. It was reported that nanofluids have good stability and 
rheological properties, dramatically higher thermal conductivities, and no penalty in pressure drop (see 
Daungthongsuk and Wongwises 13 ). A comprehensive literature on the topic of nanofluid has been discussed 
in the book by Das et al. 15 and in the papers by, Kakac and Pramuanjaroenkij 16 , Wong and De Leon 17 , Saidur 
et al. 18 , Fan and Wang 19 and very recently by Jaluria et al 20 and Mahian et al. 21 . 

Different from the above investigations, there are also many research articles on convective flow and heat 
transfer due to a shrinking sheet in recent years. Physically, there are two conditions where the flow towards a 
shrinking sheet is likely to exist, first, imposing an adequate mass suction on the boundary (Miklavcic and 
Wang 22 ) and second, consider a stagnation flow (Wang 23 ). It was reported that the analytical solutions on the 
viscous flow over a shrinking sheet with suction effects were first reported by Miklavcic and Wang 22 . It was 
observed that mass suction is necessary to maintain the flow within the boundary layer. After the earlier work by 
Miklavcic and Wang 22 , the study of flow due to a shrinking sheet was extended to other types of fluid (see Hayat 
et al. 24 and Sajid et al. 25 ). Since then, various aspect of flow and heat transfer due to a shrinking sheet has been 
investigated by researchers 26 32 . The case of unsteady flow towards a shrinking sheet was investigated by Bachok 
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et al. 33 , Rohni et al. 34 , Fang et al. 35 and Zheng et al. 36 . However, only a 
small amount of works were considered the flow and heat transfer 
characteristics due to a nonlinearly stretching or shrinking sheet 
(Rana and Bhargava 7 , Vajravelu 37 and Cortell 38 ). 

In the present study, we investigate numerically the flow and heat 
transfer over a nonlinearly shrinking sheet immersed in a nanofluid 
with suction effect at the boundary. The present study is the exten- 
sion of Rana and Bhargava 7 , to the case of shrinking sheet with 
suction effect. Numerical solutions are obtained using a shooting 
method. The effects of suction, thermophoresis, Brownian motion 
and stretching/shrinking parameters on the velocity, concentration 
and temperature profiles as well as heat transfer characteristics are 
graphically presented and discussed. 

Problem formulation. Consider a steady, laminar, incompressible 
and two dimensional boundary layer flow and heat transfer of a 
viscous nanofluid past a permeable stretching/shrinking sheet 
coinciding with the plane y = 0 and the flow being confined to y 
> 0. It is assumed that the pressure gradient and external force are 
neglected in this problem. The flow is generated by the nonlinear 
stretching/ shrinking sheet along the x-axis where x is the coordinate 
measured along the stretching/shrinking sheet. Under the boundary 
layer approximations, the governing equations for conservation of 
mass, momentum, thermal energy and nanoparticle concentration of 
this problem can be expressed as (see Rana and Bhargava 7 , Khan and 
Pop 12 , and Buongiorno 39 ) 



du dv 
dx dy 
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The boundary conditions of Eqs. (l)-(4) are 

u= XU W , v = v w , T=T W , C = C W at y = 0, 
u— >0, T—*T m , C— > Coo as _y— >co, 



(5) 



where u and v are the velocity components in the x and y directions, 
respectively, T is the temperature, C is the nanoparticle volume frac- 
tion, T w is the surface temperature, T„ is the ambient temperature, 
C w the nanoparticles volume fraction at the plate and is the 
nanoparticles volume fraction far from the plate, v w is the suction 
or injection velocity with v w < 0 for suction and v w > 0 for injection, 
t = (pc)„/(pc)f, where (pc) p is the effective heat capacity of the nano- 
particles, (pc)j is the heat capacity of the base fluid, a. = k m /(pc)f is the 
thermal diffusivity of the fluid, v is the kinematic viscosity, D B is the 
Brownian diffusion coefficient, D T is the thermophoretic diffusion 
coefficient and X is the stretching/ shrinking parameter with X > 0 for 
a stretching surface and X < 0 for a shrinking surface. The constant n 
is the nonlinearity parameter with n = 1 for the linear case and n¥= 1 
is for the nonlinear case. It is assumed that the surface is stretched or 
is shrunk with the velocity U w = ax", where a > 0 is a constant. 

Further, we seek for a similarity solution of Eqs. (l)-(4) subject to 
the boundary conditions (5) by introducing the following similarity 
transformation 7 : 



u — ax n f (fj), V- 



«"(»+!) x (n-l)/2 



n-1 
n + 1 



nf'in) 



0(n) = (T-T m )/(T w -T x ) (j>(i 1 ) = (C-C m )/(C w -C m ), (6) 



av(n+l) , ,w, 

n=y\J 2v V "-'>/ 2 

where prime denotes differentiation with respect to )}. To have sim- 
ilarity solutions of Eqs.(l)-(5), we assume 



av(n + l) oi-D/i- 
x S 



(7) 



where the constant parameter S corresponds to suction (S > 0) and 
injection (S < 0) or withdrawal of the fluid, respectively. 

Substituting Eq. (6) into Eqs. (2)-(4), where Eq. (1) is identically 
satisfied, we obtain the following ordinary differential equations: 



r + #"-^(/o 2 =o 



— 0"+fO' + NbO'<t>' + Nt(0') 2 = 0 
Pr 



The boundary conditions (5) reduce to 

/(0) = S, f'(0)=l, 0(0) = 1, 0(0) = 1 , 
/'(??)-> 0, 0(i?)-> 0, 0(i?)->O as rj-yco, 



(8) 
(9) 
(10) 

(11) 



where Pr = via is the Prandtl number and Le = v/D B is the Lewis 
number. The constant dimensionless Brownian motion parameter 
Nb and thermophoresis parameter Nt are defined as 

Nb=D B x{ - Cw ~ Cm) , Nt=D T x(T ™~ T ^ (12) 

V V To, 

The physical quantities of interest in this study are the local Nusselt 
number Nu x and the local Sherwood number Sh x which are defined 
as 



K V J w 1 co ) 



x q m 
D B (C W — Coo) 



(13) 



where k is the thermal conductibility of the nanofluid, and q w and q m 
are, respectively, the heat flux and mass flux at the surface (plate), 
given by 



<Jw = 



f d T 

Substituting (6) into (13) and (14), we obtain 



-'" [d £)„. (14) 



Re7 1/2 Mi, = - 



-^±V(0), R e -^Sh x = - ] f^f(0) (15) 



where Re x = L/^x/v is the local Reynolds number. 

It is worth mentioning that an anonymous reviewer has pointed 
out that for the linearly (i.e. n = 1) stretching/shrinking sheet, the 
exact solution for the flow field is given by 

f(tj) = S + 6(1 - e- crl ), (c = S + b>0) (16) 

with then be = X from the boundary condition/ (0) = X. Substituting 
(16) into Eq. (8) for n = 1, gives 

2 



-Sc 



X=0 



(17) 



and then 
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so that (18) gives, as it is expected, A c = — S 2 /4 < 0. 

However, for the nonlinearly (i.e. 1) stretching case (i.e. A > 0), 
there is no exact solution of Eq. (8) (see Vajravelu 37 and Cortell 38 ). 



Results and discussion 

The system of nonlinear ordinary differential equations (8)-(10) 
with the boundary conditions (11) was solved numerically using a 
shooting method. The results obtained show the influences of the 
non-dimensional governing parameters, namely suction parameter 
S, Lewis number Le, thermophoresis parameter Nt, Brownian motion 
parameter Nb and stretching/shrinking parameter A on the velocity 
profile, temperature profile, nanoparticle concentration profile, the 
local Nusselt number and the local Sherwood number. Using the 
present method, dual solutions are found by employing different 
initial guesses for the missing values of/"(0), — 0'(0) and — <j)'(0) 
where all velocity, temperature and nanoparticle concentration pro- 
files satisfy the infinity boundary conditions (11) asymptotically with 
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-3-2-1 0 1 2 
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Figure 1 | Variation of Re~ 1 l 2 Nu x with 1 for various values of S when 
n = 2, Le = 2, Nt = 0.5, Nb = 0.5 and Pr = 6.2. 




X 

Figure 2 | Variation of Re^T l ' 2 Sh x with A for various values of S when 
n = 2, Le = 2, Nt = 0.5, Nb = 0.5 and Pr = 6.2. 

two different boundary layer thicknesses. To validate the numerical 
results obtained, the present results for — 0 '(0) were compared with 
those obtained by Rana and Bhargava 7 for the case of a stretching 
surface by setting/(0) = 0 and/(0) = 1 in the boundary conditions 
(11), and neglecting the thermophoresis parameter Nt and Brownian 
motion parameter Nb, i.e. by setting Nt = 0 and Nb = 0. The present 
results were also compared with those reported by Cortell 38 by setting 
/(0) = 0 and/(0) = 1 in the boundary conditions (11) and taking 
Eckert number E c = 0 in Eq. (15) of that paper. The comparisons as 
presented in Table 1 show a favorable agreement, thus gives confid- 
ence to the results that will be reported for the shrinking case. 

The variations of the local Nusselt number Ke~ 1 ' 2 Nu x and the 
local Sherwood number Re~ 1 / 2 Sh x with stretching/shrinking para- 
meter A for S = 2.5,3 and 3.5 are shown in Figs. 1 and 2, respectively. 
From Figs. 1 and 2, we found that it is possible to obtain dual solu- 
tions of the similarity equations (8)-(10) subject to the boundary 
conditions (11). As shown in Figs. 1 and 2, it is seen that dual solu- 
tions exist for the particular suction parameter S = 2.5,3 and 3.5. 
These findings are supported by the results reported by Fang et al. 35 
for a permeable shrinking sheet that mentioning the dual solutions 
occur only with mass suction S s 2. For negative values of A, there 
exist a critical value A c with two solutions exist for A > X c , a unique 
solution obtained when A = A c and no solutions exist for A < A c . 
Based on our computations, these critical values of A c are —1.4278, 
-2.0561, -2.7986, -3.6553 and -5.7114 for S = 2.5,3,3.5,4 and 5, 
respectively. These values are presented in Table 2. 

In the following discussion, we term the first and second solutions 
in the following discussion by how they appear in Fig. 1, i.e. the first 
solution has a lower value of Re^ l l 2 Nu x than the second solution for 
a given A. For the stretching case (A > 0), dual solutions are found to 
exist for all positive values of A, to much higher values than shown in 
Figs. 1 and 2. The range A c £ A £ 2 was chosen because the solution 
shows complicated behaviors in this range. For the range A > 2, these 
quantities show monotonically increase/decrease behavior. This 
finding is in accordance with the results reported by Bachok et al. 33 
that solving the unsteady boundary-layer flow and heat transfer of a 
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Table 3 | Values of — $'(0)and — ^'(0) for different values of Sand 
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( ] second solution. 

nano fluid over a permeable stretching/shrinking sheet. From Fig. 1, 
it is seen that the values of Re x 1 / 2 Nu x which represents the heat 
transfer rate at the surface increases as the suction parameter S 
increases. This is due to the fact that increasing S is to decrease the 
thermal boundary layer thickness and in turn increase the temper- 
ature gradient at the surface. Fig. 1 also indicates that the critical 
values \k c \ for which the solution exist increase as S increases. This 
finding suggests that suction widens the region of dual solutions to 
the similarity equations (8)-(ll). 

The variation of the local Sherwood number Re.~ l l 2 Sh x with k is 
shown in Fig. 2. It is found that the first solution has a higher value of 
Re-^ 2 Sh x for a given a than the second solution. The value of 
Re x l l 2 Sh x which represents the wall mass transfer rate increases as 
S increases as illustrated in Fig. 2. The stability analysis of multiple 
solutions has been discussed in a number of studies, for example by 
Weidman et al. 40 , Paullet and Weidman 41 , Harris et al. 42 , and Ro°ca 
and Pop 43 . They have shown that the first solution is stable and 
physically relevant in practice whilst those on the second solutions 
are not. We expect that this finding holds for the present problem. 
On the other hand, the values of —O'(0) and — 0'(O) with different 
values of S for A = — 0.5 and 2 are shown in Table 3. 
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Figure 3 | Variation of Re^T l l 2 Nu x with 1 for various values of Pr when 
n = 2, Le = 2, Nt = 0.5, Nb = 0.5 and S = 2.5. 
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second solution 



012345678 

Figure 4 | Effect of the suction parameter S on the velocity profiles f (//) 
when n = 2, Le = 2, Nt = 0.5, JV6 = 0.5, Pr = 6.2 and A = -0.5 (shrinking 
surface). 

Fig. 3 shows the variations of the local Nusselt number Re^ l l 2 Nu x 
with k for several values of Prandtl number Pr and fixed values of the 
other parameters. Generally, local Nusselt number increases as Pr 
increases. From Fig. 3, it is clearly shown that for small values of Pr, 
i.e. Pr = 0.7 to 1.6, the first solution has a higher value of Re~ l l 2 Nu x 
for a given k than the second solution. However, the opposite beha- 
viors are observed for the local Nusselt number for moderate values 
of Pr, i.e. Pr = 1.8 to 6.2 as presented in Fig. 3. It is noticed that the 
second solution has a higher value of Re.~ l l 2 Nu x compared to the 
first solution for a given k. 

Figs. 4 and 5 have been plotted to demonstrate the effects of 
suction parameter S on/ (if) representing the velocity profile for both 
shrinking and stretching cases. For the first solution, it is observed 
that the velocity increases as S increases for the shrinking case as 
apparent in Fig. 4. As may be expected, it is because of the fact that 
suction cause the reduction of momentum boundary layer thickness 
and consequently enhances the flow near the solid surface. In Fig. 4, 
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second solution 
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Figure 5 | Effect of the suction parameter S on the velocity profiles f (//) 
when n = 2, Le = 2, Nt = 0.5, JVb - 0.5, Pr - 6.2 and 2 = 2 (stretching 
surface). 
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Figure 6 | Effect of the thermophoresis parameter Nt on the temperature 
profiles hjf) when S = 2.5, Nb = 0.5, Pr = 6.2, Le = 2, n = 2 and 
2 = —0.5 (shrinking surface). 

we also found that the velocity profiles have positive velocity gradient 
for the first solution, while the opposite trends are observed for the 
second solution. In Fig. 5, the fluid velocity inside the boundary layer 
decreases with an increase in S for the stretching case. This obser- 
vation occurs due to suction effect which retards the fluid motion as 
presented in Fig. 5. 

Figs. 6 and 7 give the effects of the thermophoresis parameter Nt 
on the temperature profiles 6(rf) while Figs. 8 and 9 are devoted to see 
the influences of Nt on the nanoparticle concentration profiles (f>()i). 
Figs. 6 and 7 elucidate that the temperature profiles as well as the 
boundary layer thicknesses of the thermal field increase with increas- 
ing Nt. As a result, increasing Nt is to decrease the local Nusselt 
number. Both Figs. 8 and 9 indicate that increasing Nt is to increase 
the concentration for both first and second solutions, and in con- 
sequence decrease the nanoparticle concentration gradient at the 
surface. Hence, the local Sherwood number is expected to decrease 
as Nt increases. This phenomenon is due to the thermophoresis 




012345678 
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Figure 8 | Effect of the thermophoresis parameter Nt on the nanoparticle 
concentration profiles </>{i/) when S = 2.5, Nb = 0.5, Pr = 6.2, Le = 2, 
n = 2 and 2 = —0.5 (shrinking surface). 

effect, which warm the fluid in the boundary layer. This finding is 
in good agreement with that reported in Fig. 6 in the paper by Rana 
and Bhargava 7 for the case of flow and heat transfer over a nonli- 
nearly stretching sheet in a nanofluid without suction effect. It is also 
found that the positive concentration gradient at the surface <//(0) is 
obtained for both solutions of Nt = 0.5 and 0.8 while negative con- 
centration gradient at the surface <j) ' (0) is obtained for both solutions 
of Nt = 0.1 as shown in Figs. 8 and 9, which is consistent with the 
result presented in Fig. 2. 

Figs. 10 and 1 1 are presented to observe the effect of the Brownian 
motion parameter Nb on the temperature profiles 0(1]) with the 
corresponding nanoparticle concentration profiles (f)(ii) being shown 
in Figs. 12 and 13 for both shrinking and stretching cases. These 
Figs. 10 and 11 indicate that by increasing the Brownian motion 
parameter Nb, the temperature and the thermal boundary layer 
thickness increase. This phenomenon leads to decrease the local 
Nusselt number. Different nanoparticles have different values of 
Nb and Nt. This leads to different heat transfer rate. These two 




first solution 
second solution 



Nt = 0.1, 0.5, 0.8 
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Figure 7 | Effect of the thermophoresis parameter Nt on the temperature Figure 9 | Effect of the thermophoresis parameter Nt on the nanoparticle 



profiles ffjj) when S = 2.5, Nb = 0.5, Pr = 6.2, Le = 2, n = 2 and 
2 = 2 (stretching surface). 



concentration profiles (j)(if) when S = 2.5, Nb = 0.5, Pr = 6.2, Le = 2, 
n = 2 and 2 = 2 (stretching surface). 
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Figure 10 | Effect of the Brownian motion parameter Nb on the 
temperature profiles 0(ij) when S — 2.5, Nt — 0.5, Pr — 6.2, Le = 2, n = 2 
and X = —0.5 (shrinking surface). 

parameters can be used to control the heat transfer rate in a nano- 
fluid. Also, we can see that in Figs. 10 and 11, the usual decay occurs 
to the temperature profiles for all values of Nb considered, and the 
thermal boundary layer thickness increases rapidly for large values of 
Nb. It is observed that the effect of Nb on the nanoparticle concen- 
tration profile (j>(ri) is in the opposite manner to that of temperature 
profiles 0{rf) as illustrated in Figs. 12 and 13. It is apparent from 
Figs. 12 and 13 that nanoparticle concentration is decreasing as Nb 
increasing. It seems that the Brownian motion acts to warm the fluid 
in the boundary layer and at the same time exacerbates particle 
deposition away from the fluid regime to the surface which resulting 
in a decrease of the nanoparticle concentration boundary layer thick- 
ness for both solutions. As a consequence, the concentration gradient 
at the surface increases and in turn increases the local Sherwood 
number. This finding is in accordance with the result reported in 
Fig. 4 by Rana and Bhargava 7 . The thermal conduction can be 
enhanced in two ways by the Brownian motion of nanoparticles. 




11 

Figure 1 1 | Effect of the Brownian motion parameter Nb on the 
temperature profiles 0(i/) when S — 2.5, Nt — 0.5, Pr — 6.2, Le = 2, n = 2 
and X = 2 (stretching surface). 




n 

Figure 12 | Effect of the Brownian motion parameter Nb on the 
nanoparticle concentration profiles when S = 2.5, Nt = 0.5, Pr = 6.2, 
Le = 2, n = 2 and X — —0.5 (shrinking surface). 

First, it can be enhanced by a direct effect owing to nanoparticles 
that transport heat and the second mechanism is by an indirect 
micro-convection of fluid surrounding individual nanoparticles. 
The high values of Nb means Brownian motion is strong for the small 
particle and the opposite case is applied for small values of Nb (Rana 
and Bhargava 7 ). Thus, from Figs. 10-13, it is clearly indicated that 
Brownian motion parameter provides important effect on temper- 
ature and concentration. 

Figs. 14 and 15 are depicted to examine the effects of the Lewis 
number Le on the temperature profiles 0(i]) while the corresponding 
nanoparticle concentration profiles fyiif) are presented in Figs. 16 and 
17 for both shrinking and stretching cases, respectively. It is seen 
from Figs. 14 and 15, the temperature increases as the parameter Le is 
increased, resulting in an increasing manner of the thermal boundary 
layer thickness, which consequently reduces the local Nusselt num- 
ber. However, the nanoparticle concentration profile acts in the 
opposite behavior with the increasing of Le for both solutions as 




first solution 
second solution 



0 1 2 3 4 5 

Tl 

Figure 13 | Effect of the Brownian motion parameter Nb on the 
nanoparticle concentration profiles </>(ij) when S = 2.5, Nt = 0.5, Pr = 6.2, 
Le = 2, n = 2 and X = 2 (stretching surface). 
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first solution 
second solution 



0.8 



Figure 14 | Effect of the Lewis number Le on the temperature profiles 
%) when S = 2.5, Nt = 0.5, Nb = 0.5 Pr = 6.2 and 1 = -0.5 (shrinking 
surface). 




first solution 
second solution 



Ie=l,2,5 



10 12 



14 



Figure 16 | Effect of the Lewis number Le on the nanoparticle 
concentration profiles </>{if) when S = 2.5, Nt = 0.5, JVb = 0.5, Pr 
and A = —0.5 (shrinking surface). 



= 6.2 



shown in Figs. 16 and 17. It is found that the concentration of both 
solutions decrease as Le increases. There would be a significant 
reduction in the concentration boundary layer thickness for both 
first and second solutions when Le is increased. This phenomenon 
occurs due to Lewis number effects which increases the concentra- 
tion gradient at the surface, and as a result increases the local 
Sherwood number. It is also seen that the boundary layer thicknesses 
of nanoparticle concentration for the second solutions for some 
values of Le are larger than the first solutions in both stretching 
and shrinking cases, which support the instability of the second 
solutions. The positive concentration gradient at the surface <f)'(0) 
is obtained for both solutions for small values of Le i.e. Le = 1 and 2 as 
depicted in Figs. 16 and 17. The opposite trends are displayed in the 
direction of the concentration gradient at the surface when large 
values of Le are applied for example Le = 5 as shown in Figs. 16 
and 17. It is worth mentioning that in all Figs. 4-17 presented here, 
the velocity, temperature and nanoparticle concentration profiles 



satisfy the far field boundary conditions (11) asymptotically, which 
support the validity of the numerical results obtained. 

Conclusions 

This paper presents a similarity solution of the boundary layer flow 
and heat transfer over a nonlinearly stretching/shrinking sheet 
immersed in a nanofluid with suction effect. By means of similarity 
transformation, the governing mathematical equations are reduced 
into ordinary differential equations which are then solved numer- 
ically using a shooting method. The effects of some governing para- 
meters namely suction parameter, thermophoresis parameter, 
Brownian motion parameter and stretching/shrinking parameter 
on the flow, concentration and heat transfer characteristics are 
graphically presented and discussed. The findings of the numerical 
results can be summarized as follows: 




first solution 
second solution 



0.8 



Figure 15 | Effect of the Lewis number Le on the temperature profiles 
0(ij) when S = 2.5, Nt = 0.5, Nb = 0.5, Pr = 6.2 and A = 2 (stretching 
surface). 
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first solution 
second solution 



Le = \,2, 5 



10 12 14 



11 



Figure 17 | Effect of the Lewis number Le on the nanoparticle 
concentration profiles ^(//) when S = 2.5, JVf = 0.5, Nb = 0.5, Pr = 6.2 
and 2 = 2 (stretching surface). 
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1) The local Nusselt number and the local Sherwood number 
which respectively represent the heat transfer and mass trans- 
fer rates increase with an increase in the suction parameter. 

2) The increasing of thermophoresis parameter Nt and the 
Brownian motion parameter Nb is to increase the temperature 
in the boundary layer which consequently reduces the heat 
transfer rate at the surface. 

3) A rising value in Nb and the decreasing in Nt produce a 
decrease in the nanoparticle concentration, as a result increases 
the local Sherwood number. 

4) The increase of Lewis number Le leads to an increase of the 
temperature but a decrease in the nanoparticle concentration. 

5) Dual solutions exist for a certain range of the stretching/ 
shrinking parameter a for both shrinking and stretching cases. 
For the shrinking case, the solution exists up to the critical 
values of A c (<0), while for the stretching case, the solution 
could be obtained for all positive values of X. 

6) The critical value \k c \ increases as the suction parameter S 
increases, suggest that suction widens the range of k for which 
the solution exists. 
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